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QUANTITATIVE EQUIDISTRIBUTION PROPERTIES 
OF TORAL EIGENFUNCTIONS 

HAMID HEZARI AND GABRIEL RIVIERE 


Abstract. In this note, we prove quantitative equidistribution properties for orthonor¬ 
mal bases of eigenfunctions of the Laplacian on the rational d-torus. We show that the 
rate of equidistribution of such eigenfunctions is polynomial. We also prove that equidis¬ 
tribution of eigenfunctions holds for symbols supported in balls with a radius shrinking 
at a polynomial rate. 


1. Introduction 

In [Sh74l [Z~e87 , iCdV85j . Shnirelman, Zelditch, and Colin de Verdiere proved that, on a 
compact connected Riemannian manifold (M, g ) without boundary, whose geodesic flow is 
ergodic for the Liouville measure, the eigenfunctions of the Laplacian are quantum ergodic. 
Quantum ergodicity means that, for any orthonormal basis of eigenfunctions, there exists 
a full density subsequence along which the associated microlocal lift on the unit cotangent 
bundle S*M tend weakly to the Liouville measure on S*M. The main example of ergodic 
geodesic flow is given by the geodesic flow on negatively curved manifolds. Thus, in this 
geometric context, eigenfunctions of the Laplacian are quantum ergodic. On a general 
compact Riemannian manifold ( M,g ), the geodesic flow is not ergodic for the Liouville 
measure, and the above result can be extended by using the ergodic decomposition of the 
Liouville measure — see for instance [TTTTd] . 

A natural example is the case of the rational torus T d = R^/Z^ endowed with its canonical 
metric. In this setting, the geodesic flow is not ergodic, but if one considers symbols 
dependent on x and independent of £, then an ergodic property holds for such symbols 
- see proposition O below. Using this observation, one can show that eigenfunctions 
become equidistributecQ on the configuration space [ MaRul2 , Ril3] , In this paper, we 
will give quantitative versions of this equidistribution property of toral eigenfunctions. In 
fact, these quantitative equidistribution properties are originated and motivated by some 
conjectures concerning the eigenfunctions of negatively curved manifolds. The two related 
topics we shall be concerned with are: 

• The rate of equidistribution (Theorem II.2p . 

• Small scale equidistribution (Corollary 11.51) . 


Date: March 11, 2015. 

^Note that they do not equidistribute on S*T d . 
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Also, in the last section, as an easy corollary of Zygmund’s theorem and the equidistribu- 
tion results of |MaRul2l IRil3] , we show that the quantum ergodicity holds on the 2-torus 
for L 2 symbols. 

Remark 1.1. It is worth mentioning that on the rational torus all the quantum limits of 
eigenfunctions of the Laplacian are absolutely continuous with respect to the Lebesgue 
measure [Ja971 1AM14] . In fact in dimension two, by Zygmund’s theorem, there is a 
uniform bound for the L 4 norm of all L 2 normalized eigenfunctions, and hence all quantum 
limits in this case have density functions in L 2 . This was refined by Jakobson in [[Ja97 ] 
where he proved that the density function is a trigonometric polynomial. 

Before we state our results, we fix some notations. 

Throughout the paper we denote by T d := M rf /Z rf the rational torus with the standard 
metric, and we denote by dx the normalized volume measure induced by the standard 
metric. 

Let a > 0 be some fixed positive number. For every 0 < h < 1, we consider an orthonormal 
basis {'4’f i )j=i....,N(h) °f the subspace 

Ws:=l[ 1 -afi, 1 +aft](-/» 2 A)L 2 (T d ), 

made of eigenfunctions of —ti 2 A. According to the Weyl’s law - see e.g. [ DuGu75l . one 
has N(h) ~ aA c ih 1 ^ d for some constant depending only on d. For each 1 < j < N(K), 
we denote Ej{h) € [1 — ah , 1 + ah] to be the eigenvalue corresponding to 

= mWn- 


1.1. Rate of equidistribution. Our first result states that 


Theorem 1.2. Let T d = be the rational torus with d > 2 and let a be an element 

in C°°(T rf ) (independent of h). Then, there exists some constant C a > 0 such that, for 
any orthonormal basis (V , i)j=i,...v(ft) °f l[i-aft,i+aft] ( — h 2 A)L 2 (T d ) made of eigenfunctions 


of -h 2 A, 


1 

N(h) 



al'ipjp 2 dx 



2 


< c a hl. 


The fact that this quantity converges to 0 was already observed in [MaRul2l IRil3j . The 
novelty here is that we are able to prove that this convergence holds at a polynomial rate. 
Although this is a natural question, but to our knowledge it has not been addressed in the 
literature. This result is a direct consequence of Theorem II. 41 below which is slightly more 
general. In jMaRu!2] , Marklof and Rudnick proved that equidistribution on configuration 
space also holds for eigenfunctions of a rational polygon, hence it would be natural to 
understand if one can obtain a polynomial rate of convergence in this setting. In the case 
of the torus, it would also be natural to understand what the optimal rate of convergence 
should be depending on the dimension. 
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Remark 1.3. Let us rewrite this statement using the standard convention (the non-semiclassical 
notation). We define the eigenfunctions 'ipj to be the nonzero solutions to 

—A ipj = A ji/jj, Xj > 0. 

where the eigenvalues Xj are sorted as 


0 — Ai < A2 < A3 • • • —> 00. 


Suppose (ijjj)jE n is an ONB of L 2 (T d ) made of eigenfunctions. Then, in this notation, the 
above result can be written as follows 


V(a, A) 


1 

W( A) 


E 

j:X—y/~X<.Xj <A+\/A 




0(A-3), 


where (see jDuGu75l |DiSj99| ) 

JVr(A) := )) [j : A — < Xj < X + Vxj ~ C d A^, 

with C d > 0 depending only on d. For chaotic systems, it is conjectured in the physics 
literature IFePeSGI . that V(a, A) is of order A“ 2 “. In the case of negatively curved mani¬ 
folds, the best known upper bound is 0(\ log A|" 1 ) [Ze94i ISch06 ]. We emphasize that our 
general strategy is the same as the ones in these two references: the main inputs are that 
we are able to use the semiclassical approximation for much longer times and that we have 
a better control on the error terms due to the exact formulas one has on T d . Finally, in 
the case of Hecke eigenfunctions on the modular surface, we note that the upper bound 
0(X~ 2 +e ) was proved in [LuSa95] for spectral intervals of the form [A, 2A]. 


1.2. Small-scale equidistribution. Our next result concerns equidistribution properties 
of toral eigenfunctions in balls of shrinking radius. This question is motivated by our 
recent work |HeRil4j , where we showed that on negatively curved manifolds quantum 
ergodicity holds for symbols carried on balls whose radius shrink at a logarithmic rate (see 
also |Yol31 lHan!4l ILMR15) ), and where we found some applications to L p estimates and 
the size of nodal sets. In the case of T d , one can also prove a quantitative equidistribution 
result where symbols are allowed to depend on h. This is the content of the following 
theorem which is our main result: 

Theorem 1.4. Let T d = M. d /Z d be the rational torus with d > 2. Let s > uo > 0, 
and zq > 0. Suppose a = (a^)o<s<i € C°°(X rf ) is a symbol such that for every /3 in N d , 
there exists Cp > 0 satisfying 

(1) VO < h < 1, Vx € T d , |5fofi| < Cph~ V1 ^. 

Then, there exists Hq > 0 such that, for any 0 < h < Hq, and for any orthonormal basis 
(ip{)j=i,...N(K) °f l[i-ah,i+ah}{-ft 2 A)T 2 (T*) made of eigenfunctions of-h 2 A, one has 

2 

^ 11 11 is (T d } h v ° +c*| | a/i 11 ^- s (Td) h?~ 2v °+c a h 2 ~ 2 ^ 0+L ' 1 ^, 




























4 


HAMID HEZARI AND GABRIEL RIVIERE 


where C is independent of a, uq and v\, and where C a depends only on vq, v\ and on a 
finite numbed of the constants Cp appearing in ©■ 


As was already mentioned, this result implies Theorem 11.21 by picking a independent of 
H, is i=0 and = \. This theorem also allows us to show that on the rational torus, 
the eigenfunctions equidistribute on balls whose radius shrink at a polynomial rate. More 
precisely, by choosing an, to be certain cutoff functions supported in geodesic balls of radius 
h ui , and, using an extraction argument (see for instance section 3.2 of jlleRil 1 and the 
proof of corollary 1.7 in (Han m, we get 


Corollary 1.5. Let 0 < u\ < ■ Then there exists 0 < Hq < 1/2 such that given any 

orthonormal basis {fi>n)j=i,...N{h) °f ^-[i-a.n,i+ah]{~ h 2 A)L 2 {T d ) made of eigenfunctions of 
-ti 2 A, one can find a full density subsequence A ui (h) of {1,... , N(h)} such that 


(2) VO < h < h 0 , Vx € T d , Vj € A U1 (h) : 


. Sb(x,H- 1 ) \'K( x )\ 2dx . 
ai - Vol(B(x, hL 1 )) -° 2 


where the constants 01,02 > 0 are independent of h, x, and j a,d B(x,h Ul ) denotes the 
geodesic ball of radius h Vl centered at x. 


The fact that (. A Ul (H))o<n<h 0 is °f full density exactly means that 


lim 


IA u 1 {h)\ 

N{h) 


= 1 . 


In other words, this statement says that there exists a large proportion of eigenfunctions 
where the average value of the square of eigenfunctions in shrinking balls of radius h vi , 
are uniformly bounded by two constants. We note that the corollary provides a subse¬ 
quence of density 1 subsets that works uniformly for every point on the torus. If we had 
considered only one fixed point xq in T d , we would have obtained a critical exponent of 
size instead of the exponent appearing here. The exponent appearing 

in this statement is probably not optimal and it is plausible that this exponent can be 
improved by using for instance methods like the ones used by Bourgain in |Bol3] , Here, 
our proof relies only on tools from ergodic theory and semiclassical analysis. 


Remark 1.6. Using this corollary and the strategy of (HeRil4j . we can in fact improve 
Sogge’s L p estimates for toral eigenfunctions jSo88j . However the L p bounds we obtain 
using this method are not better than the upper bounds proved in |Zy74[ Bo93 . IBol.31 
lBoDe!4] . 


1.3. Quantum ergodicity for L 2 observables. Motivated by the recent question raised 
by Zelditch on L°° quantum ergodicity |Zel3| . we mention the following nice consequence 
of the quantum ergodicity property on the 2-torus: 

Corollary 1.7. Let T 2 = R 2 /Z 2 be the rational 2-torus. Then, for any orthonormal 
basis °f ^[\-ah,i+ah}{~h 2 A)L 2 (T 2 ) made of eigenfunctions of—h 2 A, there 

2 The /3 involved in the constant depends on the choice of uq and v\. 
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exists a full density subsequence A(h) of {1,2,..., N(h)} such that, for all a(x) £ L 2 (T 2 ) 
(independent ofh), 


( 3 ) 


lim 

h-^o,j£A(h) 



a(x)\ip{(x)\ 2 dx 



a(x)dx. 


The important point in this statement is that convergence holds for any observables in 
L 2 (T 2 ), and not only C° ones (as it is usually the case in quantum ergodicity statements). 
In particular, it holds for the characteristic function of any measurable subset of T 2 . In 
fact, Zelditch conjectured in }Zel3j that (|3|) holds for any a in L°°(M ) provided (M, g ) is a 
negatively curved manifold. As will be explained in section [U this corollary follows directly 
from a classical result of Zygmund |Zy74| combined to the above quantum ergodicity 
property. We emphasize that we do not need all the strength of the above theorems and 
that this result could be in fact deduced directly from the results in }MaRul2l !Ril3] . 

We will now give the proof of Theorem ll.4l from which all the other results follow. As in the 
case of negatively curved manifolds, we will first prove a result on the rate of convergence 
of Birkhoff averages. Then, we will implement this result in the classical proof of quantum 
ergodicity, and we will have to optimize the size of the different remainders to get our 
results. 

Remark 1.8. After communicating this note to Zeev Rudnick, he informed us that part 
of the above results can in fact be improved using methods of more arithmetic nature. 
The proof presented here only makes use of standard tools of Fourier analysis, and it is 
modeled on arguments similar to the ones used to prove rate of quantum ergodicity on 
negatively curved manifolds | Ze94 . Sch061 . 


2. Convergence of Birkhoff averages 


We start with the following proposition which gives us the rate of equidistribution for 
observables depending only on the x variable: 


Proposition 2.1. Let T d = R. d /Z d with d > 2. There exists C d > 0 such that, for every 
a in C°°(T d ), one has 


I S*T d 


1 f T f 

— / a(x + tf)dt — / a(x)dx 

T Jo J T d 


dxdf < 


C d \\a 


|2 

lL 2 (T d ) 


Proof. Let a be a smooth function on T d . We write its Fourier decomposition a := 
Y^kez d a k£k, where e k (x ) := e 2mk ' x . We set 


V[a,T) 


j j — I a(x + tf)dt — j a(x)dx 
h d Js d_1 T J 0 Jjd 


2 

dfdx. 
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First, we perform integration in the x variable and we find that 


V(a,T) = -j ]T \a k | 


k£Z d -{ 0} 


r T ' 2 

e 2M^ dt 


/ 0 




Now we would like to estimate the integral in each term of the above sum. Using the 
spherical symmetry first, and then calculating the dt integral, we get 


rT i 2 

e 2inkt ^dt 


to 




3 2i7r||fc||^i^ 


di 


sin 2 (7 tT||A:||6) 




Js d -1 ^Pll 2 ^ 

By using the spherical coordinates and putting £i = cos<^ , 0 < y? < n, we obtain 


f sin 2 (7rTj|fcj|^i) 
/gd-i vr 2 ||A:|| 2 ^ 


d£ = C 


f'* sin 2 (7rT||fc|| cosy?) 
o vr 2 ||fc|| 2 cos 2 y? 


(sin ip) d 2 dip. 


where the constant C is the value of the integral with respect to the remaining spherical 
variables. The change of variable s = cos y>, turns this last integral into 


J0 TT 2 \\k\\-S 2 


We then split this integral into integrals over [0,5] and its complement [5,1], where 0 < 
5 < 1. Clearly 


[ 1 <_ - 

Is vr 2 ||fc|| 2 s 2 (V j ~7 t 2 P| 


p|| 2 5 2 ' 

To estimate the integral on [0, 5], we use the substitution u = 7tT||A;||s. Hence, since d > 2 
we get 


r s sin 2 (7rT||fc|ls) 
to vr 2 ||A:|| 2 s 2 


(V 7 ! - s 2 ) d ~ 3 ds < 


r-7rT||fc||«5 gin 2 n 


< 


00 sin 2 u 


7r||A;||\/l — 5 2 do 
T 


7r||A:| | \/l — 5 2 Jo u2 
T 

2||A:||VT^‘ 


du 


u* 


du 


Therefore, by choosing 6 = -^== and T > 2, we get 


a 




2 r'lh/ii 2 

a k \ ^ ° l|a|| L 2 (T d ) 




for some uniform constant C > 0. 


□ 
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3. Proof of Theorem 11.41 


We fix a in C°°(T rf ) that potentially depends on h , even if we omit the index h in order to 
alleviate the notations. We also suppose that a belongs to a nice class of symbols. More 
precisely, there exists v\ > 0 such that, for every a in N rf , one can find C a > 0 such that 

(4) V.x € T d , |0“a| < C a h~ U1 ^. 

Without loss of generality, we will also suppose that a is real valued. We set a := a — 
f Td adx. Our goal is to give an upper bound on the following quantity: 

i m r r 

v h, 2 (a) = ——^2 / a(x)\il) 3 h (x)\ 2 dx- / adx 

Jl d J T d 


3.1. Applying Egorov’s theorem. We rewrite the previous expression as follows 


1 7V(K) 

i4 ’ 2(o> =7v(ft) x; 

3 = 1 


4* 


One of the main differences with the negatively curved case treated in |Ze941 ISch06L lHan!4l 
lHeR,il4| is that we can consider much longer semiclassical times using the fact that we 
are on T4 Precisely, we fix uq > 0 and T = T(K) := h~ u °. We introduce the averaged 
operator 




_ ithA _ ithA 

e 2 ae 2 dt. 


Using the fact that 'ipl is an eigenmode for every 1 < j < N(ti), one can write that 


N(h) 


Vh ’ 2{a) = Wh)^ \{^A(T,n)H) L2 

V ' 3 = 1 


Using the Cauchy-Schwarz inequality, we find that 

1 JV(ft) 


V K ,_{a) < — (^ h ,A[T,h)H\ 


L 2 


Recall now that (V ; /i)j=^ is an orthonormal basis of the space 

^ = l[i-^, M] (-^A)L 2 (T d ). 

We will now take advantage of the fact that we are on the torus. Precisely, we have 


( 5 ) 


14 , 2 ( 0 ) < 


N(h) 


E 


\\A(T,h)e k || 


L 2 


: (27rh \| k | |) 2 E [1— ah, l-\-ah] 


where e k (x) := e 2mk - x . We write the following exact formula 


(A(T, h)e k ) (x) = 


rT(h) 


T(h) 


E 


(XriG 


p£Z d — {0} 


'e p (x + 2Trhtk)dtj e k (x), 
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where a := a p e p (x). We now fix s > We note that, uniformly for x in T d and 

t in R, one has 


ith{2 / K) ||p|| 

a p e 2 e p (x + 2nhtk) — a(x + 2n kht) 

pdlA— { 0 } 


< 2n 2 ht'^^ ||p|| 2 |d p | < c s ht\\a\\H s , 

p^O 


for some constant c s > 0 depending only on cL and s. Implementing this in our upper 
bound <[5|), we obtain that 
( 6 ) 


Vs ,2 (a) < 


N(h) 


E 


1 


/ T d 


T(h) 


k G : ( 2nh \ \ k 11) 2 E [ 1— ah, l+ah] 

where the constant in the remainder is independent of a. 


rT[h) \ 

J a(x + 2irkht)dt J dx+\\a\\ 2 HS 0(h 2 2v °), 


3.2. Trace asymptotics. In order to compute the previous expression, we proceed as 
in [ DuGu75j - see |DiSj99| (Ch. 11) for a semiclassical version. We will in fact follow 
the presentation of Prop. 1 in [Sch06] and we will take advantage of the fact that we are 
working on T4 


Regarding Q, we now have to estimate 


14,2 (a) 


1 

N[h) 


E 

k€jj d : (2nh\\k \ |) 2 E [1— oth, l-\-ah] 



bfr(x, 2irkh)dx. 


where we set 

1 f fi “" 0 

J 0 a 

with 0 < \i < 1 a smooth cutoff function which is equal to 1 in a small neighborhood of 
1 and is 0 outside a slightly bigger neighborhood, say outside [1/4,4], We fix a smooth 
function p > 0 in the Schwartz class 5(M), which is > 1 on the interval [—a,a]. We also 
suppose that p has compact support, say that the support is included in [—1/8,1/8]. 


(x + t£)dt 


hM := x i (11 Cl 1 2 ) 


Remark 3.1. In order to construct such a function, one can start from a nonzero smooth 
even function / > 0 which is compactly supported in [—1/16,1/16] and take p to be the 
inverse Fourier transform of A^f * / with Aq > 0 large enough. We note that the function 
p satisfies p'( 0) = 0. 


We can then write that 


v hAa) = hWlEi) Wl , 2M)ix , 


N(h) 


fee z d 




i4TT 2 T\\k\\ 2 h 2 


r)e ft e ft 


/ bn(x, 2nkh)dxdr. 
Jf d 


and thus 
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Thanks to the Poisson summation formula, we get 


(7) V H , 2 (a) = 


1 


N(Jl) JY d 


E 

dez d 


( 27 vhy 


■ [ p(t) 

■J RxR d 


Tdisik-u-g-i , , 

e ft o/j(x,£; )drd^ ax. 


We will now make use of the stationary (and non-stationary) phase lemma. To do so, 
we fix / in and we denote by <pi(r,£) to be the phase function of the above oscillatory 
integral. We observe that, for l / 0, one has ||d^^|| > |j/|| — 2r||£|| > ||2|| — 1/2, for r in 
the support of p and ||£|| 2 in the support of xi- For l / 0, we introduce the operator 

h d^ipi-dg 


P, := - 


i || d^ipi 


12 ' 


We perform N integration by parts using this operator and we find that, for every x in 
T d , 


1 


(27 rh) 1 


[ P(r)t 
J RxR d 


1 -ft- b h (x,0dTd£ 


< C 




-1/2) 


N ’ 


for some uniform constant C > 0 that depends only on p , a and d. Using the upper 
bound 0 and taking N large enough in the previous equation, we get 


(8) Vn ,2 {a) < 


1 


N(h) 


'T d 


(27T H) d XxR d ^ ^ 


e i Am n 1) b h (x,C)dTd e ) dx + O(h 2 ^- U0 -^). 


We now disintegrate the measure d£ along the energy layers |||£|| 2 — 1 = 77} (or in other 
words we use the coarea formula), to write for every x in T d , 

r+oo 


[ P(r) 

J RxR d 


(iifir-D , , 
e" ft b h (x, 


r r+oo tE 

£,)drd£ = / / p(r)e l ~ (bn) {x, E)dEdr , 

Jr J -1 


where 


(h)(x,E):= f b h (x,£)dL E (£). 

J\\t\\ 2 -l=E 


We can now use the stationary phase formula and the fact that p'(0) = 0 - see for 

instance |Zwl2| (Ch. 3). Precisely, we find that 

(9) 


/ 

J RxR d 


r(liqr-l) 


p{r)e l ft bft(x, ^)drd^ = p(0)2irh 


(iien 2 =i 


b h (x,0dL 0 (0 + O(h 2 ^-^)) . 


As before we take N large enough (depending only on vq , and d) in the stationary phase 
lemma to ensure that the remainder is of order O(h 2 ^ 1 ~ uo ~ u i)). We give emphasis that the 
remainder term is of the form ^ 2 ( 1 _!/ o-^i) anc [ no t 2(i'o+*'i). this is due to the fact that 
p'( 0) = 0 and to the symmetry of the phase function - e.g. Theorem 3.17 of [Zw 12]- 

Now by the Weyl’s law, we know that N(K) ~ aCdbJ~ d , for some constant depending only 
on d - see for example Ch. 11 of |DiSj99| . Thus, combining ©, © and Q, we have that 

rh~" o 


V h , 2 (a)<C 0 f [ [ a(x + t£)dt 

Jf d a Jo 


d£dx + C'||a||^ft 2 " 2l/ ° + O{h 2{1 - U0 ~ Vl) ), 


2 
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where Co depends only on d and on the choice of p, and where C' s depends only on s. 
We also note that the constant in the remainder depends only on finitely many of the C a 
appearing in Q. 


3.3. The conclusion. We can now apply Proposition 12.11 and we finally find that 
Vn, 2 (a) < CoC d \\a\\ 2 L2{Td) h uo + 0(h 2 ^~^) + C' 0 \\d\\ 2 HS{Jd) h 2 - 2u °, 
which concludes the proof of Theorem 11.41 


4. Proof of Corollary 11.71 

The proof of this result is a direct consequence of the quantum ergodicity property on the 
2-torus, of Zygmund’s theorem on the L 4 norms of the eigenfunctions on the 2-torus |Zy74| , 
and of the Banach-Alaoglu theorem. 

Let (V’fi)j=i,-..A r (?i) b e an orthonormal basis of ^-[i- a h,\+ah] (—^i 2 A)L 2 (T 2 ) made of eigen¬ 
functions of -h 2 A on the rational 2-torus. Then by Theorem 11.21 there exists a full 
density subsequence of A (h) of {1,2,..., N(h)} such that for all a(x) € C°(T 2 ) 

(10) lim [ a(x)\^i(x)j 2 dx = [ a{x)dx. 

h-^o,j£A(h) Jj2 Jj2 

We refer to section 15.4 in [Zwl2j for the details of the extraction argument. We want to 
show that (11011 holds for all a € L 2 (T 2 ). To prove this, we first note that by Zygmund’s 
theorem wm , there exists a uniform constant A such that 

f |V4( X )1 4 dx < A 4 . 

J T 2 

Thus, the sequence T := {\ipf i \ 2 )i<j<N(h),o<h<i is bounded in L 2 . By the Banach-Alaoglu 
theorem, it is relatively compact for the weak-* topology on L 2 (T 2 ). On the other hand, 
by (1101) . T has at most one weak-* limit in L 2 , and that is the constant function 1. This 
proves the corollary. 
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